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ABSTRACT 


A  series  of  reacting  flow  tests  of  a  solid  fuel  ramjet 
were  conducted  with  different  air  inlet  swirl.  Cold  flow 
measurements  of  the  air  flow  at  the  exit  plane  of  the  tuhe- 
in-hole  injector  were  used  to  determine  the  Swirl  Number  for 
each  configuration.  Regression  rates  of  an  HTPB  type  fuel  at 
low  air  mass  flux  were  found  to  increase  with  swirl  vane 
angle  (and  Swirl  Number)  up  to  30  degrees. 
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INTRODUCTION 


A.  SOLID  FUEL  RAMJET 

Ramjet  technology  has  taken  three  distinct  approaches  to 
meet  the  requirements  of  an  ever-expanding  mission  responsi¬ 
bility  'placed  upon  weapon  systems  of  today.  These  three 
concepts  are  the  Liquid  Fuel  Ramjet  (LFRJ),  Solid  Fuel  Ram¬ 
jet  (SFRJ)  and  the  Ducted  Rocket  (DR).  Each  system  has 
distinct  advantages  over  the  other.  The  LFRJ  has  the 
capability  to  operate  over  widely  different  levels  of  thrust 
and  efficiency  through  use  of  variable  fuel  control.  The 
SFRJ  is  lighter,  more  compact,  and  is  much  less  complex 
than  the  LFRJ  due  to  the  absence  of  fuel  tanks,  fuel 
delivery  and  control  mechanisms.  The  DR  retains  some  of 
the  benefits  of  the  SFRJ,  but  must  carry  oxidizing  materials 
within  the  fuel  to  maintain  combustion. 

Because  of  the  simplicity  with  which  the  SFRJ  operates, 
and  its  relative  compactness,  it  is  desireable  to  develop 
a  simple  thrust  control  for  the  solid  fuel  ramjet  motor.  To 
understand  how  this  problem  may  be  approached,  it  is  first 
necessary  to  have  a  workinq  knowledge  of  the  basics  of  ramjet 
technology. 

As  stated  earlier,  the  ramjet  is  a  very  simple  engine. 
Since  there  is  no  compressor,  the  ramjet  uses  supersonic 
velocities  to  compress  the  air  flow.  Figure  1  depicts  the 
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Figure  1.  Schematic  of  Simple  SFRJ 
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basic  SFRJ.  The  three  components  of  the  basic  SFRJ  are  the 
supersonic  diffuser  or  inlet,  the  combustor,  and  the 
exhaust  nozzle. 

As  supersonic  flows  enter  the  SFRJ  diffuser,  several 
significant  events  take  place.  Standing  shocks  are  present 
at  the  air  inlet.  An  increase  in  static  temperature  and 
decrease  in  flow  velocity  occurs  through  the  shock  waves. 
Combustion  of  the  solid  fuel  increases  the  mass  flow  slightly 
and  the  temperature  is  again  increased.  This  larger  mass 
flow  at  higher  temperature  is  then  expanded  through  the 
exhaust  nozzle.  A  detailed  description  is  given  by  United 
Technologies  Corporation,  Chemical  Systems  Division  TRef.  1]. 
It  is  the  increase  (or  decrease)  of  the  fuel  prolysis  rate  or 
regression  rate  (r)  that  determines  the  thrust  and  effects 
the  efficiency  of  the  SFRJ. 

There  are  many  methods  that  can  be  used  to  affect  the 
regression  rate  of  a  solid  fuel.  One  is  to  change  the  fuel 
composition.  Another  is  to  change  the  environmental  effects 
on  the  fuel,  such  as  the  chamber  pressure  <Pc),  inlet  air 
temperature,  and  the  mass  flux  through  the  port  of  the  fuel 
grain.  A  typical  regression  rate  formula  is  given  by  Boaz 
and  Netzer  [Ref.  2] 

t  =  c  Pc0'51  G° ’ u 1  T°*33  (1) 

for  motors  not  using  bypass  air  where 


c  =  constant 

Pc  =  chamber  pressure 

G  =  mass  flow  per  unit  area  in  the  fuel  port 
Ti  =  inlet  air  temperature 

The  SFRJ  is  sel  f-thrott- 1  ing  through  the  effect  of  G  on 
r.  However,  the  self-throttling  is  not  sufficient  for  many 
missions.  Alternate  means  for  varying  r  are  therefore 
desired.  Changing  G  can  be  accomplished  by  either  changing 
m  (such  as  in  bypass  designs)  or  by  changing  the  flow  area. 
Variable  bypass  is  one  possible  means  for  controlling  r*. 

Ko  ("Ref.  3]  investigated  the  effects  of  inlet  air  swirl 
on  the  regression  rate  of  the  solid  fuel  grain.  Swirl 
imparts  an  angular  component  of  mass  flux  into  the  flow  field 
which  can  effect  the  residence  time  of  the  flow  by  increasing 
the  effective  length  of  the  combusion  section  of  the  ramjet 
motor.  It  can  also  affect  the  convective  heat  transfer  to 
the  fuel  surface  and  the  mixing  rates  of  fuel  and  oxygen.  It 
was  found  that  swirl  may  be  an  effective  means  for  regression 
rate  control,  although  only  limited  data  were  obtained. 

B.  SWIRL  NUMBER 

The  introduction  of  swirl  into  a  flow  causes  difficulty 
in  characterization  of  the  flow  field.  There  is  no  longer 
only  the  axial  component  of  the  flow  but  now  radial  and 
tangential  components  as  well.  The  Swirl  Number  (S)  is  one 
way  to  characterize  the  flow.  Lilley  and  others  TRef.  4? 
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have  done  considerable  work  with  swirling  flows  and  have  set 
forth  the  non-dimensional ized  swirl  number  (for  specific 
assumptions  about  the  inlet  swirl  behavior)  as 


S  = 


where 


<5 

0 

Gx 

d/2 


Gx  (d/2) 

=  tangential  flux  of  swirl  momentum 
=  axial  flux  of  momentum 
=  equivalent  inlet  radius 


This  thesis  investigated  the  effects  of  varying  the  inlet 
air  flow  angle  on  the  combustion  process  of  a  non-bypass 
SFRJ.  This  was  accomplished  by  using  tube-in-hole  injectors 
with  vanes  of  0,  15,  30,  and  45  degrees.  Cold  flow 

measurements  were  made  to  determine  inlet  flow  angles,  and 
reacting  flow  tests  were  conducted  to  determine  the  effects 
on  fuel  regression  rate. 


II.  DESCRIPTION  OF  APPARATUS 

A.  RAMJET  MOTOR 

The  ramjet  motor  used  in  this  investigation  was  identical 

. 

with  one  used  m  earlier  inves lqations  at  he  Naval  t 

I 

Postgratduate  School  fRefs.  2,3].  Figure  2  is  a  schematic  of  i 

the  head-end  assembly,  step  insert  section,  fuel  qrain  and  1 

nozzle. 

The  head-end  assembly  used  a  wedge  to  turn  the  primary 
inlet  air  from  90  degree  plenum  dumps.  The  igniter  torch  and  j 

I 

and  ignition  gas  were  injected  into  the  recirculation  zone  | 

just  aft  of  the  step  insert. 

The  step  insert  section  acted  as  the  air  inlet,  which 
was  modified  with  a  tube-in-hole  injector.  The  step  height  j 

I 

allowed  for  sudden  expansion  of  the  flow  into  the  combustion 
section  as  depicted  in  Figure  3.  This  sudden  expansion 
allows  a  flame  to  be  stabilize  behind  the  step  by  creating  a 
recirculation  zone.  The  swirl  elements  and  tube  were 
fabricated  from  stainless  steel.  Each  swirl  blade  was  a 
cambered  airfoil  attached  at  the  hub  with  silver  solder.  The 
vane  tips  were  ground  for  minimum  clearance  at  the  tube  inner 
wall,  and  then  silver  soldered  along  the  downstream  side  of 
the  blades.  Figure  4  shows  the  tubes  with  the  swirlers  used 
in  this  experiment.  Tubes  with  15,  30,  and  45  deg.  (from 


15 


UNSTEADY  REATTACHMENT 
ZONE 


Figure  4.  15,  30,  and  45  Degree  Swirl 

Blades  Mounted  Inside  Tubes 
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the  axial  centerline)  swirl  vane  elements  and  a  tube  without 
a  swirler  were  used. 

CSD-18818  (an  HTPB  based  system)  was  the  solid  fuel  used 
for  the  hot  runs.  The  average  length  and  internal  diameter 
of  the  fuel  grains  used  were  12.0  and  1.7  inches 
respectively.  The  fuel  grain  was  placed  between  the  head 
assembly  and  an  aft  mixing  chamber.  The  latter  contained  the 
exhaust  nozzle  and  was  connected  with  steel  rods  which  bolted 
to  the  head-end  assembly.  The  entire  motor  was  mounted  on 
the  thrust  stand.  Figure  5  is  a  photograph  of  the  thrust 
stand  a  the  test  cell.  The  control  panel  was  used  to 
initiate  the  runs  and  the  Hewlett-Packard  9836  Computer  and 
3054A  Data  Acquisition/Control  system  sequenced  events  during 
the  runs. 
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Figure  5.  Experimental  Thrust  Stand  and  SFRJ 
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Figure  6.  Schematic  of  Air  and  Gas  Supply  System 
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III.  EXPERIMENTAL  PROCEDURES 


A.  CALIBRATION 

Calibration  of  pressure  transducers  was  conducted 
at  the  beginning  of  each  test  day  utilizing  a  dead-weight 
tester.  Table  1  is  a  list  of  the  transducers  calibrated  and 
the  calibration  pressures/load. 


TABLE  1 

TRANSDUCER  CALIBRATION  PRESSURES/LOAD 


Pa 

500 

psi 

Pc 

100 

psi 

Ph 

100 

psi 

Phf 

400 

psi 

Pho 

500 

psi 

Pif 

250 

ps  i 

F 

100 

lbf 

B.  DATA  '  XTRACTION 

A  Honeywell  1508B  Visicorder  was  used  to  obtain  analog 
readings  from  pressure  and  thrust  transducers.  Chromel- 
alumel  thermocouples  were  used  for  measurements  of  all  gas 
temperatures.  A  Hewlett-Packard  9836  Computer  and  3054A 
Automatic  Data  Acquisition/Control  Stytem  were  utilized  to 
extract  digital  information. 


C.  HOT  FIRING 


Fuel  grain  length  (L),  average  port  diameter  (Dp)  ,  and 
initial  weight  (W)  were  measured  for  each  run.  The  inlet  was 
then  mounted  on  the  ramjet  motor  and  preparations  were  made 
for  the  ensuing  runs. 

Air  and  ignition  gas  pressures  were  calculated  using  the 
choked  "flow  formula  (assuming  Cd  =  0.97) 


Pressures  were  set  using  dome  loaders  and  mass  flow  rates 
were  set  using  the  computer  system  as  the  controlling  device. 
The  computer  was  then  readied  for  data  acquisiton  and  each 
test  was  initiated  from  the  control  panel.  Table  2  lists  the 
inlet  geometry  and  test  condintions  for  each  run. 

Upon  completion  of  each  run,  the  fuel  grain  was  weighed 
and  weight  loss  recorded.  Actual  motor  burn  time  (tb),  air 
time  (ta),  ignition  time  (ti),  pre-ignition  chamber  pres¬ 
sure  and  thrust  (Pa  and  Fa),  and  average  chamber  pressure 
and  thrust  (Pc  and  F)  were  obtained  from  the  traces  recorded 
on  the  Visicorder.  The  change  in  fuel  grain  weicjht  is  be 
expressed  as 

iW  =  J  <d|  -  d‘)Lpf  (4) 

Solving  for  the  average  final  port  diameter, 


»in/-n-vnn^iTrv»T  1-1  •  mn  T  70  OWTnT  ATT  TA  X TkTT'TI 


TABLE  2 


INLET  GEOMETRY  AND  TEST  CONDITIONS 


Run 

Tube  Length 

Swirl  Angle 

ta, ti, tb,  i 

(  in.  ) 

(deg. ) 

( sec.  ) * 

1 

2.5 

0 

6,1,  6,4 

2 

2.5 

0 

6,1,  6,4 

3 

2.5 

15 

6,1,  6,4 

4  * 

2.5 

15 

6,1,  6,4 

5 

2.5 

30 

6,1,  6,4 

6 

2.5 

30 

6,1,  6,4 

7 

2.5 

45 

6,1,  6,4 

8 

2.5 

45 

6,1,  6,4 

Remarks : 


Unsuccessful  first  two  attempts.  Cooled  motor, 
reweighed  grain,  cleaned  ignition  qas  choke  and 
torch.  Third  try  successful. 

pre-ignition  airflow  time,  ignition  time,  burn  time  and 
purge  time  respectively 


'■iiVdi'Ayt 


and  substituting  into 


(6) 


yeilds  the  average  fuel  regression  rate. 

Average  values  for  m(air),  m(02),  m(C  H4  ) ,  ,  and 

m(ign  fuel)  were  then  calculated  from  the  digital  printout 
for  the  run.  The  fuel  mass  flow  rate  was  calculated  by 


m 


f 


A  Wf 

"V 


(7) 


The  Naval  Weapons  Center  (NWC)  China  Lake,  Ca., 

Propellant  Evaluation  Program  (PEPCODE)  was  used  to  obtain 
theoretical  values  for  the  adiabatic  combustion  temperature 
and  combustion  gas  properties,  Y  and  R,  using  Pc,  Ti, 
heats  of  formation,  and  elemental  composition  by  weight  of 
the  fuel  grain  and  injected  gases  used  during  each  run.  Ko 
[Ref.  3]  gives  a  complete  explanation  of  the  procedures  used 
to  obtain  the  values  needed  for  data  reduction. 

D.  COLD  FLOW  PROCEDURES 

A  hemispherical,  three-hole  stagnation  pressure  probe 
was  used  to  measure  flow  velocities  at  the  exit  plane  of  the 
tube-in-hole  injector.  The  probe  was  connected  by  tubing  to 
a  pressure  transducer.  The  transducer  in  turn,  was  connected 
through  a  controller  to  a  voltmeter  where  stagnation  pres¬ 
sure  readings  in  millivolts  were  displayed.  A  Scanivalve 
allowed  for  cycling  through  the  three  pressure  readings  to 
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find  the  angle  at  which  the  flow  was  eminating  from  the 
tube-in-hole  injector  exit  plane.  This  was  accomplished  by 
finding  the  maximum  value  for  the  center  port  and  nulling  the 
left  and  right  values.  The  velocities  between  the  tube  and 
the  inlet  wall  were  also  measured. 

The  probe  could  be  rotated  about  its  own  axis.  It  was 
mounted  on  a  vernier  translation  stand  to  allow  for  movement 
across  the  flow  in  0.05  in.  increments.  The  probe  itself 
could  be  translated  on  its  own  axis  as  well  as  rotated. 

This  facilitated  taking  measurements  in  the  y-z  plane  with 
rotation  about  the  y-axis.  Figure  7  is  a  photograph  of  the 
probe  mounted  on  the  thrust  stand  in  this  configuration. 

With  the  angle  of  flow  taken  at  each  grid  point,  the  probe 
was  then  rotated  from  the  y-axis  to  the  z-axis.  Figure 
8  depicts  the  probe  head  and  the  degrees  of  freedom  in  which 
it  could  be  moved. 

The  pressure  readings  in  millivolts  and  the  angles  at 
which  the  the  maximum  values  occurred  were  then  input  into  a 
program  that  caluclated  x,  y,  and  z  components  of  the  flow 
velocity,  and  transformed  them  from  the  Cartesian  coordinate 
system  to  the  cylindrical  coordinate  system.  This  produced 
axial  and  tangential  velocity  components  (assuming  that  the 
net  radial  component  was  egual  to  zero).  The  computer 
program  and  examples  of  values  calculated  are  given  in 
Appendix  A. 
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DEGREES  OF  FREEDOM 
Y-Z  TRANSLATION 
Z-AXIS  ROTATION 


VIEW  FROM  Z-AXIS 


Y-Z  TRANSLATION 
Y-AXIS  ROTATION 


Figure  8.  Pressure  Probe  Head  and  Degrees  of  Freedom 


Upon  completion  of  characterization  of  the  flow  at  each 
point  in  the  grid,  the  amount  of  swirl  in  the  flow  was  then 


characterized  using  Swirl  Number  (S)  as  previously  stated 
in  equation  2.  The  tangential  flux  of  momentum  is 


_  it  p  u  w  (d/2) 
—  — 2 —  o  mo 


and  the*  axial  flux  of  momentum  is 


_  _  TT  p  u  (d/2j 2  (1  -  (g/2) 2 ) 

Gx  -  —  ° 

Making  the  same  assumptions  as  Lilley  TRef.  4]  of  solid  body 
rotation  (a  flat  upstream  axial  velocity  profile,  and  an 
increasing  swirl  velocity  from  0  at  the  hub  to  at  the 

outer  edge),  the  analysis  yeilds  the  Swirl  Number  given  as 


(1  -  (g/2) 2) 


(10) 


where 


is  the  ratio  of  the  maximum  tangential  velocity  to  axial 
velocity  at  the  exit  plane  of  the  tube-in-hole  injector. 

In  an  attempt  to  more  accurately  characterize  the  swirl 
in  the  present  apparatus,  a  mass  averaged  momentum  ratio  was 
calculated  by  taking  the  mass  average  tangential  velocity 


I  m  w 
I  m 
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IV. 


RESULTS  AND  DISCUSSION 


A.  INITIAL  TESTS 

A  series  of  initial  hot  runs  was  completed  to  test  the 
integrity  of  the  tube-in-hole  injector  design.  During  this 
phase  i't  was  found  that  initially  designed  tubes,  which 
extended  0.25  inch  beyond  the  aft  end  of  the  inlet,  did  not 
allow  ignition  of  the  motor.  It  was  suspected  that  swirling 
flow  was  interferring  with  the  recirculation  zone,  destroy- 
Ong  the  flame  stabilization  area.  The  tubes  were  lengthened 
so  that  the  swirl  dump  was  1.5  inches  beyond  the  inlet  step, 
which  solved  the  problem.  This  length  was  considerably  less 
than  the  reattachment  length  of  approximately  3.0  inches. 

Burn  time  for  the  motor  was  initially  set  for  10  seconds, 
which  almost  resulted  in  for  burn-through  to  the  wall.  Burn 
time  was  therefore  reduced  to  8  seconds,  which  worked  well 
for  non-swirl  test  runs.  Burn-through  occurred  at  8  secs, 
with  30  degrees  of  swirl,  so  the  burn  time  was  further  re¬ 
duced  to  6  seconds  and  the  remaining  hot  runs  were  conducted. 

B.  COLD  FLOW  ANALYSIS 

A  summary  of  results  of  the  cold  flow  studies  is  given  in 
Table  3.  There  was  an  increase  in  the  difference  between  the 
maximum  and  mass  averaged  axial  velocities  from  8%  for  no 
swirl,  to  43%  for  45  degrees  of  swirl.  The  same  was  true 
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for  the  difference  between  the  maximum  and  mass  averaged 
tangential  velocities,  with  21%  difference  at  15  degrees  and 
43%  difference  at  45  degrees.  This  could  have  been  a 
function  of  the  increased  blockage  within  the  tube.  The 
velocities  between  the  tube  and  inlet  wall  increased  as  the 
vane  angle  increased,  which  was  an  indication  that  blockaqe 
was  occurring.  The  Swirl  number,  as  calculated  by  Equation 
10  (based  on  solid  body  rotation),  increased  with  vane  angle 
to  30  degrees  and  then  decreased.  The  mass  averaged  value 
(Equation  12  divided  by  Equation  13)  increased  continually 
with  vane  angle. 

TABLE  3 

COLD  FLOW  RESULTS 


Blade  angle 
(deg) 

0 

15 

30 

45 

u ( max ) 

( ft/ sec) 

475 

616 

759 

820 

w(max ) 

( ft/sec) 

0 

309 

521 

525 

u ( avg ) 

( ft/sec) 

434 

509 

560 

462 

u( outer  flow) 

( ft/sec) 

384 

393 

503 

589 

w(avg ) 

( ft/sec) 

0 

243 

345 

296 

S(avg) 

0.0 

0.478 

0.616 

0.640 

S 


0.0 


0.335 


0.649 


0.544 


C.  FUEL  GRAIN  PHYSICAL  CHARACTERISTICS 


The  fuel  grains  used  without  swirl  exhibited  smooth 
burning  surfaces  with  the  port  diameter  slowly  increasing  to 
a  point  approximately  2.0  inches  beyond  the  step.  There¬ 
after,  the  port  diameter  was  approximately  constant.  The 
grains  used  for  the  15  degree  swirl  showed  the  same  charact¬ 
eristics.  However,  four  areas  that  corresponded  with  the 
four  vanes  of  the  swirl  element  exhibited  slightly  indented 
burn  regions  extending  from  2  inches  beyond  the  step  to  4 
inches,  slowly  blending  into  a  concentric  surface.  The  30 
degree  grains  showed  a  marked  increase  in  the  indented  burn 
regions,  and  they  extended  to  5  inches  beyond  the  step.  The 
45  degree  grains  were  even  more  pronounced  in  the  indented 
burned  region  but  not  as  elongated.  The  angles  of  the 
indented  burned  regions  corresponded  to  the  respective  swirl 
vane  angles.  Table  4  lists  the  physical  characteristics  of 
fuel  grains  pertaining  to  the  tests. 

D.  COMBUSTION  PROPERTIES 

Measured  guantities  for  each  hot  run  are  listed  in  Table 
5.  A  plot  of  the  regression  rate  vs  swirl  vane  angle  is 
given  in  Figure  9  and  vs  S  and  S(avg)  in  Figure  10. 
Regression  rate  increased  rapidly  for  small  amounts  (15 
degree  vane  angle)  of  swirl.  Further  increases  in  swirl  vane 
angles  (or  Swirl  Number)  had  much  less  effect  on  r.  The 
Swirl  Number  (Eguation  10)  correliated  reasonably  well  with 


regression  rate.  S(avg)  appeared  to  be  in  auestion  for  the 
45  degree  vane,  perhaps  due  to  the  large  amount  of  blockage. 

A  typical  value  of  G  (mass  flux)  in  an  SFRJ  is  0.6 
lbm/in2-sec  without  bypass.  Fifty  percent  bypass  would  then 
yeild  a  fuel  port  G  of  0.3.  With  r ^  G*4,  an  approximately 
25%  reduction  in  r  could  be  achieved.  A  30  degree  swirl,  in 
the  present  configuration,  would  yeild  only  about  half  of 
this  change  in  r.  However,  it  does  appear  that  variable 
swirl  can  be  an  effectiv  method  of  regression  rate  control  in 
the  SFRJ. 

TABLE  4 

FUEL  GRAIN  PHYSICAL  CHARACTERISTICS 


Run  # 

Vane 

Angle 

L 

Dp 

Di 

Dth 

Ath 

AW 

(deg) 

(  in) 

(  in) 

(  in) 

(  in ) 

(sq.  in) 

( gm) 

1 

0 

12.010 

1.750 

0.92 

0.943 

0.6483 

197 

2 

0 

12.000 

1.749 

0.92 

0.943 

0.6604 

200 

3 

15 

12.008 

1.747 

0.92 

0.943 

0.6554 

219 

4 

15 

12.014 

1.759 

0.92 

0.943 

0.6735 

222 

5 

30 

12.010 

1.754 

0.92 

0.943 

0.6597 

233 

6 

30 

12.060 

1.740 

0.92 

0.943 

0.6525 

230 

7 

45 

11.929 

1.750 

0.92 

0.943 

0.5721 

221 

a 

45 

11.915 

1.746 

0.92 

0.943 

0.6648 

220 

TABLE  5 


MEASURED  COMBUSTION  PROPERTIES 


Run# 

ti 

( sec) 

tb 

( sec ) 

ta 

( sec ) 

Pa 

(psia ) 

r 

( in/sec ) 

Pc 

(psia) 

Ti 
(o  R) 

1 

r  0.82 

7.27 

6.0 

50.3 

0.0236 

100.2 

1103 

2 

0.84 

7.26 

6.1 

49.9 

0.0240 

107.8 

1125 

3 

0.86 

7.19 

6.0 

49.5 

0.0267 

105.7 

1125 

4 

0.78 

7.20 

6.1 

50.3 

0.0266 

107.8 

1119 

5 

0.73 

7.21 

6.0 

51.1 

0.0277 

108.6 

1122 

6 

0.80 

7.22 

6.1 

51.1 

0.0274 

109.3 

1119 

7 

0.79 

7.18 

6.8 

50.7 

0.0263 

107.8 

1122 

8 

0.87 

7.27 

6.0 

49.1 

0.0276 

103.8 

1130 

m 

air 

( lbm/ sec) 

m 

fuel 

m 

°2 

m 

CH 

4 

m 

if 

m 

tot 

1 

0.503 

0.0598 

0.0190 

0.0046 

0.0033 

0.593 

2 

0.506 

0.0608 

0.0190 

0.0047 

0.0335 

0.593 

3 

0.495 

0.0672 

0.0191 

0.0047 

0.0034 

0.589 

4 

0.510 

0.0680 

0.0193 

0.0048 

0.0034 

0.604 

5 

0.508 

0.0713 

0.0191 

0.0047 

0.0034 

0.606 

6 

0.513 

0.0702 

0.0194 

0.0049 

0.0033 

0.611 

7 

0.503 

0.0602 

0.0194 

0.0048 

0.0034 

0.591 

8 

0.487 

0.0667 

0.0192 

0.0049 

0.0036 

0.581 
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REGRESSION  RATE  VS  VANE  ANGLE 
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REGRESSION  RATE 

0.02300.0235  0.0240  0.0245  0.0250  0.0255  0.0260  0.0265  0.0270  0.0275  0.0280 


REGRESSION  RATE  VS  SWIRL  NUMBER 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


It  has  been  found  that  the  regression  rate  increases 
significantly  with  small  amounts  of  swirl  in  the  flow  of 
inlet  air,  but  that  larger  amounts  of  swirl  yield  less 
effect.  The  Swirl  Number  S  (as  given  by  Eguation  10),  at 
the  low  mass  flux  values  of  approximately  0.25  which  were 
used  in  this  investigation,  correlated  reasonably  well  with 
regression  rate. 

It  is  recommended  that  the  effects  of  Swirl  Number  on 
regression  rate  be  investigated  with  increasing  values  of 
mass  flux  to  1.0  lbm/in  -sec  and  with  metalized  fuels. 
Improvement  of  the  45  degree  swirl  element  design  should  be 
also  attempted,  as  well  as  an  investigation  of  the  proper¬ 
ties  of  combustion  with  swirl  anoles  between  30  and  45 
degrees  and  beyond. 


I 
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AEEEHQI*  A 

PROGRAM  FOR  COLO  FLOW  ANALYSIS  OF  SWIRLING  INLET  AIRFLOWS 

SWIRL  NUMBER  PROGRAM 
LT.  WILLIAM  H.  CAMPEELL,  JR. 

NAVAL  POSTGRADUATE  SCHOOL,  MONTEREY 

THIS  PROGRAM  IS  USED  TO  COMPUTE  THE  AVERAGE  SwIRL  NUMBER  OF  A 
SWIRLING  FLOW  OF  A  TUBE  IN  HOLE  INJECTOR  FOR  A  SOLID  ROCKET 
RAMJET  AND  TRANSFORM  THE  CARTESIAN  VELOCITY  COMPONENTS  INTO 
CYLINDRICAL  VELOCITY  COMPONENTS. 


DEFINITIONS 


Tim 

= 

INLET  AIR  TEMPERATURE 

SONVEL  (I) 

= 

SONIC  VELOCITY 

DENS  (I  ) 

= 

AIR  DENSITY 

Yll.KI 

= 

Y-CARTESIAN  COORDINATE  OF  THE 

EXIT 

PLANE 

Z( J,1 > 

= 

Z-CARTESIAN  COORDINATE  OF  THE 

EXIT 

PLANE 

K  1  {  J  »  K  ) 

RADIUS  TO  POINT  Y(N),Z(M) 

Y  (  I,  J,K) 

= 

PRESSURE  VECTOR  IN  MILLIVOLTS 

ABOUT 

THE  Y-AXIS 

R  Z  (  I,  J,K) 

= 

PRESSURE  VECTOR  IN  MILLIVOLTS 

ABOUT 

THE  Z-AXIS 

ALFA!  I  ,  J  ,K  ) 

= 

ANGLE  MEASURED  OF  R Y ( N , M ) ,  POSITIVE 

ABOVE  THE 

£ZTA<  I  ,J,K) 

= 

ANGLE  MEASURED  OF  RZ(N,M),  POSITIVE 

LEFT  OF  Y 

All  I, J  ,K) 

= 

X  COMPONENT  GF  R Y ( N »  M  ) 

a2  (  I ,  J , K ) 

= 

X  COMPONENT  OF  RZ(N,M) 

9 


B(  I » J  »  K) 

C  (  I » J  *  K) 

xco  ypt  i  *  j  »  k)  = 

RTCU  I  «JiK)  = 

RXZ(IfJtK)  = 

HXV(I,J,K)  = 

oAMMA  (  I,  J  ,M  = 
DELTAl XVJ ,K)= 
PHKIfJtKl  = 
TC 0.>1P  (  I,  J  ,K)  = 

RC  OMP I  I, J  «K)  = 
?TX(I,J»K)  = 
P  T  T  { I  f  J«  K  )  = 

J  (  I  *  J  t  K) 
w  (  I  *  J  •  K)  = 

jUMl 


SUM2 

SJM3 


Y  COMPONENT  OF  RZ(N,M) 

Z  COMPONENT  OF  RYIN.M) 

TOTAL  X  COMPONENT  OF  PRESSURE  IN  MILLIVOLTS 
TCTAL  VECTOR  OF  PRESSURE  IN  MILLIVOLTS 
IN  MILLIVOLTS 

TOTAL  VECTOR  COMPONENT  OF  PRESSURE  IN  THE  X-Z 
PLANE  IN  MILLIVOLTS 

TOTAL  VECTOR  COMPONENT  OF  PRESSURE  IN  THE  X-Y 
PLANE  IN  MILLI VOLTS 
ANGLE  BETWEEN  Y=0  ANO  RYZ 
ANGLE  BETWEEN  Y=0  AND  Ri 

ANGLE  3ETWEEN  X  =  3  ANO  RXY,  SWIRL  ANGLE  OF  FLOW 
TANGENTIAL  COMPONENT  OF  STAGNATION  PRESSURE 
IN  MILLIVOLTS 

RADIAL  COMPONENT  OF  STAGNATION  PRESSURE  IN  MILLIVOLT 

X -COMP ON ENT  OF  STAGNATION  PRESSURE  IN  PSI 

TANGENTIAL  COMPONENT  OF  STAGNATION  PRESSURE  IN  PSI 

AXIAL  VELOCITY  COMPONENT 

TANGENTIAL  VELOCITY  COMPONENT 

SUMA  T I  ON  OF  THE  GRID  REPRE SEN  TAT  I VE  AREA 

MULTIPLIED  By  THE  GRID  POINT  TANGENTIAL 

VELOCITY  W I  1 1  J  t  K ) 

SUMA  T I  ON  OF  THE  GRID  REPRESENTATIVE  AREAS 
SUMATICN  OF  THE  GRID  REPRESENTATIVE  AREA  MULT. 

BY  THE  GRID  POINT  AXIAL  VELOCITY  J(I,J,K) 
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JM AX( I  i 


WMAXl I J 

AM  (  I  I 

All  (  I ) 

S  A  V  ( I  ) 
SMAX(  I  ) 
G  (  I  ) 


THE  MAXIMUM  AXIAL  VELOCITY  FOR  THAT  RESPECTIVE 
SWIRL  VANE  ANGLE 

THE  MAXIMUM  TANGENTIAL  VELOCITY  FOR  THAT 
RESPECTIVE  SWIRL  VANE  ANGLE 
THE  AVERAGE  TANGENTIAL  VELOCITY  FOR  THAT 
RESPECTIVE  SWIRL  VANE  ANGLE 

THE  AVERAGE  AXIAL  VELOCITY  FOR  THAT  RESPECTIVE 

SWIRL  VANE  ANGLE 

THE  RATIO  OF  THE  AW(I)  TO  AU(I) 

THE  SWIRL  NUMBER  ACCORDING  TO  LILLEY  ET  AL 
THE  RATIO  OF  THE  WMAX (  I)  TO  UMAX!  I) 


DELCARATI  ON  OF  MATRICES  FOR  RESPECTIVE  COMPONENTS  OF  THE  PROGRAM 

REAL  Y (  1 » 8 ) »  Z ( 8 1 1 ) »  TII4I,  THETA ( 4 )  »  SONVEL l 4 ) »  DENS ( 4 ) , AW < 4)  , 

♦  RY ( 4*  8  ? 8 )  «  ALFA! 4,8*8)*  RZ44,8,8),  B ETA ( 4 ,8 , 8 ) , AU ( 4 ) , 

+  XCOMP ( 4»  8 , 8 ) «  Ri  ( 8  t  8 ) *  RYZC4,8,8)»  RTOT < 4 , 8 , 8 ) , S UM1 1  4 ) » 

♦  PHU4.8.B).  RCCMP(4,8t8),  PTX(4,8,3>,  UC4,8,8),  SUM214) 

REAL  PTT (4 «  8*  8 )  »  W(4,3,8),  G(4),  AREA<8,3),  SUM3I4),  SUM4I4), 

♦  AI  (  4t  8  1 8  )  *  A2 ( 4»  8 »  8 ) »  B(4,d,3>,  C ( 4 , 8 » 8 ) t RXZ ( 4, 8 , 8 ) ♦ 

+  R  X Y ( 4  *  8 » 8 )  «  DELT  A ( 4  f  8  » 8 )  »  GAMMA! 4  » 8  »  8  ) »  TCGMP < 4 , 8 , 8 ) , 

+  MM  AX!  4 ) »  UMAX (4) «  SAV!4J,  SMAX!4) 

INCREMENTATION  OF  THE  Y  AND  Z  AXIS 
DO  l  N  = 1 . 8 
R  EAD ( 5 »  10)  Y ( I *  N ) 
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Z (N»  1)*Y  (l,N) 

CONT  INUE 
00  81  J  =  1, 8 

R£AO(  5  t  400  )  (  ARE  A (  J  t  K  )  t  K.  =  1 , 8  J 
CONTINUE 

COMPUTATION  OF  SONIC  VELOC ITY, OEMS  ITY,  ANO  INPUT  OF  VANE  ANGLE 
DO  2  1*1,4 

READC5  t  20)  TI  (  I)  ,  THE  TA  ( I  ) 

SONVEL (  I>  =  SGRT(1  .4*32 .2*53 . 3*T I { I )  ) 

OENS(I )=14.7/(53.3*TI(I ) ) 

CONT  INUE 

THIS  PORTION  OF  THE  PROGRAM  INPUTS  ANO  OUTPUTS  DATA  RECORDED 
DURING  THE  COLD  FLOW  ANALYSIS  OF  DIFFERENT  SWIRLER  VANE 
GEOMETRIES. TH  VALUE  (8,8)  FOR  EACH  ANGLE  IS  THE  COMPONENT 
MEASURED  BETWEEN  THE  INLET  WALL  AND  THE  TUBE-IN-HOLE  INJECTOR 
OUTER  WALL. 

DO  4  1=1,4 
SUM  1 (I) =0 
SUM2 (I) =0 
SUM3 ( I ) =0 
SUM4 ( I ) =0 
UMAX ( I ) =0 


WMAX  (I ) =0 
DO  5  K  =  1 1 3 
DO  6  J  =  l»  8 

READ(5,30)  RY(I,J,K)»  AL FA < I  * J, * >  ,  RZ(I,J,K),  8ETA{I,JfK) 
AL(  I » J,K)=RYC I* J,K)*CDS (ALFA( I , J,K) *3. 141/18  0) 

A2( IfJ#K)=kZ( I»J»K)*CDS(8ETA{I»J»K) *3. 141/180) 

B ( I  *  J*K )  =  RZ  (  I. J,K)*SIN <BETA(I  , J, K ) *3 . 141 / 1 80 ) 

Cl  I  *  J%K»=RY ( I t J, K ) *S IN ( ALFA ( I. J,K 1*3.141/ 180) 

XCOMP!  I »J*K)  =  A1(I,J,K)  +  A2< I » J « K ) 

Rl( J,K)=SQRT( YClfKI**2+Zt J,l>**2) 

RYZ (IfJ,K)=SGRT(C!I,J,K)**2+6< I,J,K >**2) 

RXZ  ! I,J,K)=SQRT(C( I( J,K)**2+XCQMPI I,J,K)**2) 

RXY  (  I»J»  K)  =SQRT  (XCOMP! It J»K)**2+B(I  ,JtK)**2) 

RTQT  (It  JtKl  =SCRT!  XCOMP ! I , J , K ) * *2+6!  I , J , K ) **2+C 1 1 , J , K ) **2 ) 
IF(  (Y(l tKJ.EQ.O)  .AND.!Z(J»l).EQ.O))  THEN 
OE  LTA( I , J  *  K J =0 

ELSE  IF  ( (Z( Jtl) .GE.0.05)  .AND.  (Y(  I  »K)  .tO.O  )  )  THEN 

DELTA! I, J ,K) =90 

ELSE 

OELTA( I ,J  »K) =(ATAN(Z(J*I)/Y(1 » K) ) ) *1 80/ 3 • 14 
END  IF 

IF  (  e  ( I  *  J.KJ.EQ.OJ  THEN 
GAMMA!  1 1 J  *  K )  =0 
ELSE 

GAMMA! I , J ,K) =( ATANtC! I , J , K ) / 8 ( I , J , K ) ) ) *  1 8 0/ 3. 14 
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yvV-V-V 

-  *  j  »•».  • 


END  IF 

IF  (XCOMPI I,J,K).EQ.O>  THEN 
PHI(I«JtK)>0 
ELSE 

PH  1(1, J,K)  =  ( ATAN(B< I, J,K)/XCOMP( I, J,K) ) >*180/3.14 
ENO  IF 

TCOMPU,  J,K)=RYZ(I,  J,K)*CDS(  (90- (DELTA!  I  ,  J  ,  K )  ♦GAMMA!  I  ,  J  ,K  )  )  ) 

♦  *3.141/180) 

RCQM  P(I»J»K)=RYZ(I»J»K)*SIN((90-!DELTA{  I , J , K ) +GAMMA < I , J ,K )  )  ) 

+  *3.  141/180) 

PTX (  I,J,K)=XCOMP( I , J ,K )/.0 155 

UII , J,K)=SQRT( 2402. 7*TI ( I )* ( 5*( I !PTX ( I , J , K) / 14. 7) ♦ l ) ** .2857 1 
♦  l)  1  ) 

PTT  (  I,J,K)=TCCMP(I  ,J,K  J/.0155 

W(I  , J,K)=SQRT(2432.7*TI ( I ) * ( 5* Cl ( PTT ( I , J ,K ) / 14.7 )  +  1 ) **  .285  71 

♦-1J  )  ) 

SUM  1(1)*  SUNK  I)+(AREA(J»K.)*W(I  »  J»K)  ) 

SUM  2(1) =S  JM2 (  I ) +AREA ( J , K ) 

SUM  3 ( I ) =SUM3 ( I )♦( AREA( J ,K) *U(  I  ,J,K) ) 

IF  (U(IfJ,K).GT.UMAX(I ) )  THEN 
UMAX!  I)  3U  (  I,  J,  K) 

ELSE 

UM  AX ( I )=UMAX (I ) 

ENO  IF 

IF  (W(IfJfK).GT.WMAX(I)}  THEN 


it  j.kj 

ELSE 

WM  AX(  I )  =*WMAX  ( I ) 

ENO  IF 
CONTINUE 
CONTINUE 

AW(  I)=SUM1(  I)/SUM2(I) 

AU(  n*-SUN3(  I ) /SUM2 ( I  I 
S AV ( I ) =AW( I )/ AU(  I) 

S(  I  )  =  W  MAX ( I )/ UMAX ( I ) 

SMAXU  )*(G(l)/2l/(l-3(I  1**21 


CGNTINJE 

WRI  T  E  (  6  »  41 0  ) 

(UMAX!  I) 

t  1*1,4) 

WRITE(6 ,411  ) 

( WMA  X ( I ) 

•> 

r-« 

II 

mRITE(6 .412) 

(  AW (  1 1  , 

1  =  1,4) 

*RI T  E(6  » 413 1 

(  AU(  I), 

1  =  1,4) 

WRI T E (6  » 414 ) 

(SAV(  I), 

1*1,4) 

WKI  TE( 6 . 415 ) 

(  SMAX(  I) 

.  1*1,4) 

If 

M 

QO 

o 

o 

«RITE( 6,401 

WRI TE( 6,501 

THE TA(  I  ) 

( Y( 1  ,L), 

L* l , 3  ) 

WRITE( 6,55) 

0G  9  J=l,3 

WRITE (6,60) 

Z( J  ,L)  , 

(U(I,4,K),K=l,3) 

CONTINUE 
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K 

U~ 

m 


i 


3  CONTINUE 

00  12  1*1,4 
WR I TE4 6*70)  THETA(I) 

WR  I  TE4  6 ,50 )  4  Y4 1 ,L),L=1,8) 

WR  I  TE4  6 ,55) 

00  11  J-1,8 

WRITE (6*60)  Z ( J  *1 ) *  (H(I.J*KI»K-1.8) 

11  CONTINUE 

12  CONTINUE 

DO  14  1*1,4 

WR  ITE  (  6,80  )  THETA(I) 

WR  I  TE  (  6,50  )  (  Y4  1,L)  *L-1,8) 

WR ITE ( 6,55 ) 

00  13  J*1 , 8 

WRITE (6, 90)  Z(J,1),  (RY ( I  , J  ,K)  ,  K*l,8) 

13  CONTINUE 

14  CONTINUE 

00  16  1*1,4 
WR  I  TE(  6,80)  THETA(I) 

WR  ITE4  6,50)  (  Y  ( 1  ,U,l*l,3  ) 

WRITE! 6,55) 

00  15  J*1 , 8 

WRITE (6,120)  Z(J,i),  4PHI(I,J,K),  K*l,8) 

15  CONTINUE 

16  CONTINUE 
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17 

13 


21 

22 


23 

24 


00  18  1*1,4 

WRITE (6,100)  THETA(I) 

WRITE! 6,50)  ( Y(  1«L) ,LS 1 « 8 ) 

WRITE (6,55) 

00  17  J=1 , 8 

WRIT  E  (6, 90  )  Z(J,1>,  (  RZ  (  L  *  J  ,  K)  ,  K=l,8) 
CONTINUE 
CONTINUE 
DO  22  I  =  1,4 

WR ITE (6,110)  THETA(I) 

WR  ITE  (6,50)  €  Y I  1,U  ,L*1,8) 

WRITE (6,55) 

00  2  1  J=  1 , 8 

WRIT  E (6, 120  )  Z(J,1),  ( AL  FA ( I,J,K),  K=l,8) 
CONTINUE 
CONTINUE 
DO  24  1*1,4 

WRITE (6,130)  THETA(I) 

WR ITE ( 6,50 )  ( Y( l,L),  L*l,8) 

WRITE (6,55) 

DO  23  J=l , 8 

W  RI T  t { 6, 120  )  Z(J,1),  <  SET  A( I , J , K) ,  K=i,3) 
CONTINUE 
CONT  INUE 
DC  32  1  =  1,4 
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WRITE !6,140)  THETA(l) 

WRITE  16(50)  4  Y { 1,0,  1*1, 8i 
WRITE (6,551 

00  31  J*  L  *  8 

WRITE(6,150)  Z!J,1),  (C(I,J,K),  K*  1 , 8  ) 

31  CONTINUE 

32  CONTINUE 

DO  34"!  *1,4 

WRITE ( 6*160)  THETA(I) 

WRITE  (6,50)  |Y(ltLI,  L=l,8) 

WR ITE ( 6,55 ) 

DO  33  J* 1 , 8 

WRIT  E (6, 1 50 )  Z  ( J , 1 ) ,  ( XCOMP ( I  ,  J , K) ,  K*l,8) 

33  CONTINUE 

34  CONTINUE 

00  36  1*1,4 

WR ITE (6,170)  THETA(I) 

WRITE  (  6,50)  (Y(  1,0, 1=1, 8) 

WRITE  <6,55) 

DO  35  J  =  1 » 8 

WRI T  E (6, 1 50 )  Z(J,l),  ( RTOT ( I , J , K) ,  K=I,8) 

35  CONTINUE 

36  CONTINUE 

DO  38  1=1,4 

WRITE! 6, 180 1  THETA! I ) 
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WRITE  (6*50)  IV(UL),  1  =  1,3) 

WR ITE ( 6,55 ) 

00  37  J  s  1 , 8 

WRITE (6, 120 )  Z(J ,1),  (PHI(I»J»K) 
CONTINUE 
CONT  INUE 
00  42  1=1,4 

WR ITE l 6, 190)  THETA(I) 

WHITE (6,50)  ( Y( 1, L) ,  L  =  1 , 8  ) 

WRITE (6,55) 

00  41  J=l »  8 

WRITE (6, 150  )  Z(J,1),  (Al( I • J, K)  , 
CONTINUE 
CONT  INUE 
00  44  1=1,4 

WR ITE ( 6,200)  THETA(I) 

WR  ITE  (  6,50 )  <Y(l,l),  L=l,8t 
WRIT  E (6, 55 ) 

DO  43  4=1,8 

WRITE(6,150)  Z(J,l),  (  A2 ( I , J , K)  , 
CONTINUE 
CONT  INJE 
00  46  1=1,4 
WRITE( 6,210)  THETA(I) 

WR  I  TE(  6,50  )  <  Y  (  1  »L  )  ,  L=l,8) 


WRITE ( 6,55 ) 

00  45  J*1 y 8 

WRIT  E (6, 150 )  ZlJfll.  !B(I,J,K),  K*l,3) 

45  CONTINUE 

46  CONTINUE 

00  48  1=1,4 

WRITE  16,220)  THETA(I) 

WR  ft  E  (  6, 50 )  f  Y  (  1,  L)  ,  L  =  l,8) 

WRITE!  6,55 ) 

DO  47  J  =  1 ,8 

WRITEI6, 150 )  Z(J,1),  !RyZ(I,J,K),  K=l,8) 

47  CCNTINUE 

48  CONTINUE 

DO  52  1  =  1,4 

WR ITE (6,230)  THETA! I ) 

WRITE ( 6,50 )  ( Y( 1, L) ,  L=l,8) 

WRITE (6,55) 

00  51  J=l , 8 

WR ITE (6,150)  Z!J,1),  (  RX  Z ( I , J, K)  ,  K=l,8) 

51  CONTINUE 

52  CONTINUE 

OG  54  '=1,4 

WRITE  (6,240)  THETAU) 

WRITE (6, 50)  ( Y ( l,L) ,  L=l,8) 

WRITE!  6, 55) 
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00  53  J»l,8 

WRITE (6, 150  )  Z(J»1),  (RXY(I,J,K),  K=l,8) 

53  CONTINUE 

54  CONTINUE 

DO  57  1=1,4 

WR ITE ( 6,250 )  THETA(I) 

WRITE (6,50)  ( Y ( 1 , L J ,  L=I,8) 

WR  HE!  6,55  ) 

DO  56  J=1 , 8 

WRI TE  (6, 120  )  Z(J,i),  (DELTA! I,J ,K) ,  K=l,8) 

56  CONTINUE 

57  CONTINUE 

DO  59  I  =1,4 

WRITE  (6,260)  TrlETA(I) 

WRITE!  6*50)  (Y(  1,L),  L  =  l,8J 
rfRITE (6,55) 

DO  58  J* 1 , 8 

WRI T  E (6, 120 )  Z(J,1),  <  GAMMA! I , J ,K) ,  K=l,8) 

58  CONTINUE 

59  CONTINUE 
DO  62  1=1,4 

WRITE!  6, 270)  THETAU  ) 

WRITE (6,50)  ( Y ( 1,L),  L=l,8) 

WRITE (6,55) 

DO  6 l  J= 1 , 8 
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WRIT E (6,150)  Z(J.l),  4TCOMP( I,J ,K) ,  K=l,3) 
CONTINUE 
CONTINUE 
00  64  1  =  1,4 

WRITE! 6,280)  THETA(I) 

WR ITE  4  6,50 )  ( Y  4 i,L) ,  L  =  1 ,8) 

WRITE (6,55) 

0  0V6  3  J  =  1 ,8 

WRIT  E (6, 150 )  Z(J,1),  (RC0MP4 I, J ,K) ,  K=L,3) 
CONTINUE 
CONTINUE 
DO  66  1=1,4 

WRITE!  6,290 )  THETA(I) 

WRITE (6,50)  4  Y ( 1 , L) ,  L  =  1 , 8  ) 

WRITE (6,55) 

00  65  J=1 , 8 

WRIT  E  (6,60)  Z(0,l),  (PTX(l,4,KI,  K.=  l,8) 
CONTINUE 
CONT INUE 
00  68  1=1,4 

WR I TE ( 6 » 3 10 )  THETA(I) 

WRITE(6,50)  (Y(l,L),  L  =  l,8) 

WRITE ( 6,55) 

00  67  J  =  l » 8 

WRIT  E  (  6,  601  Z(J,i),  (PTT(I,J,K),  K=l,3) 
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67  CONTINUE 

68  CONTINUE 

WRITE (6,330) 

MR ITE (6,50)  ( Y( I, L),  L=i ,8) 

WR ITE ( 6, 55 ) 

00  73  J=I,8 

WRIT  E (6,60)  Z 1 J, 1 J •  (RlfJtKJ,  K=l,8) 


73 

CONTINUE 

10 

FORMAT!  F6.3) 

20 

FORM  AT ( 2F6.3) 

30 

FORMAT!  4F10. 3) 

40 

FORMAT( •I,,'U(I»J»K)  FOR  SWIRL  VANE 

BLADE 

ANGLE' ,F5.1, 

>.'/  ) 

50 

FORMAT!  9X,'Y*  ,8F3.3) 

55 

FORMAT!  ' O' »  3X •  *Z'  ) 

60 

FORMAT!  F6.3,4X,8F  8.2) 

70 

FORMAT!  •  1*  ,  'M(  I  »J,K)  FOR  S*IRL  VANE 

BLADE 

ANGLE* ,F5. I, 

♦  .’/) 

80 

FORMAT (  '  I','RY(  I  *  J, K)  FOR  SWIRL  VANE 

i  BLADE 

ANGLE' ,F5.I 

+'.'/) 

90 

FORMAT!  F6.3,4X,8F  8.3) 

IOC 

FORM  AT (  '  I', '  RZ  1  I,J,K)  FOR  SWIRLVANE 

BLADE 

ANGLE' , F  5* I , 

♦  .  '/  ) 


113  FORMAT  ( *  l ' » ' ALFAl  I, J, K )  FOR  SWIRL  VANE  BLADE  A  NGLE  * , F  5 


120  FORM AT I F6*3,4X»8F8«0) 

130  FORMAT !*l*»'BETAl  I,  J, K )  FOR  SWIRL  VANE  BLAOE  ANGLE*,  F5. 
♦I,* . '/) 

140  FORMAT ! *i*»'C(  I , J  «K  )  FOR  SWIRL  VAN  E  BLADE  ANGLE* ,F5. 1 , ' 
♦  .*/) 

150  FORMAT ! F6«3»4X  »  8F  8.3 ) 

160  FORMAT l*l*»'XCOMPCI»J,K)  FOR  SWIRL  VANE  BLADE  ANGLE*, F5 

*  * 

+  .1, '  .*/  J 

170  FORM  ATI *1*,'RT0T( I, J, K)  FOR  SWIRL  VANE  BLADE  ANGLE*, F5. 
♦  1,'  . */) 

180  FORMAT ! 'l*,*PHl(I,J,K)  FOR  SWIRL  VANE  BLADE  ANGLE', F5.1 

190  FORMAT! ' 1*, 'Alt  I, J,K)  FOR  SWIRL  VANE  BLADE  ANGLE',  F5. 1, 
«■'.*/) 

200  FORMAT!  •  1*,  •  A2(  I,  J,K)  FOR  SWIRL  VANE  BLADE  ANGLE*, F5.1, 
♦  *.'/) 

210  FORMAT! « l*, 'B! I, JtKI  FOR  SWIRL  VANE  BLADE  ANGLE  * , F 5. L , • 

♦  .'/ I 

220  FORMAT!  •  l*, *RYZ! I  ,J,K)  FORSWIRL  VANE  BLADE  ANGLE*, F5.1, 

♦  *.*/) 

230  FORMAT! • 1* , ' P X2 ! I , J,K)  FOR  SWIRL  VANE  BLADE  ANGLE' ,F5.1 

♦, ». */) 

240  FORMAT! • l« ,*RXY( I  ,J»K)  FOR  SWIRL  VANE  BLADE  ANGLE', F5.1 
♦  ,  *.  • /) 

250  FORMAT!  •  1* , • DELTA ! I , J , K)  FOR  SWIRL  VANE  BLADE  ANGLE', F5 


+  .1,  ' .'/ ) 

260  FORM  AT(  '  1*  , '  GAMMA  <  I  ,  J  ,  K )  FOR  SWIRL  VANE  BLADE  ANGLE*  »  F5 

♦  .l.  '  .'/  ) 

270  FORMAT! • 1* , • TCOMP (I , J, K)  FQRSWIRL  VANE  BLAOE  ANGLE*, F5. 

♦  1,*.  '/I 

280  FORMAT!  •  1  • , • RCOMP  (  I  , J , K J  FORSwIRL  VANE  BLADE  ANGLE', F5. 

♦1,* . */) 

290  FORMATt • l* , *PTX! I  ,J,K)  FOR  SWIRL  VANE  BLADE  ANGLE*  ,F5 • I 

♦  ,  *.  * /) 

310  FORMAT! • i*,'PTT< I ,J,K)  FOR  5'WIRL  VANE  BLADE  ANGLE*, F5.1 

♦  ,'.*/) 

STOP 

330  FORMAT!  *  l*, 'R1!J,K)  FOR  ALL  EXPERIMENTS'/) 

400  FORMAT!  8F8. 4) 

410  FORMAT!  ' 1* , •  UMAX(I):  •  , 4{ F6. 2 , • , • , 2X ) ) 

411  FORMAT!//, '  WMAXII):  *  ,  4 ( F6 . 2 , * , • , 2X )) 

412  FORMAT!//,*  AW(I):  • , 4<F6.2, • , ' , 2X1) 

413  FORMAT!//,'  AU ( I) :  * , 4( F6 . 2 , ' , * , 2X) ) 

414  FORMAT!//, *  SAV(I):  •  ,  4  !  F6. 3  ,  •  ,  •  ,2XH 

415  FORMAT!//,*  SMAX(I):  •  ,4  l  F6 . 3  ,  '  ,  •  ,  2X  ») 

END 

SENTRY 

'JMAX(I):  475.00,  616.81,  759.39,  820.37, 

WMAX(I):  0.  00,  309.46,  521.23,  525.  85, 
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AW1I1 


0.00*  243.41*  345.32*  296.17* 


AU(I):  434.54*  509.01*  560.96*  462.46* 


SAV( I ) :  0.000,  0. 

478,  0 

1.616* 

0.640 

» 

SMAXm:  0.000,  0 

.335, 

0.649* 

0.544, 

u(  i*  j*kV 

'  FOR  SWIRL  VANE 

i  BLADE 

ANGLE 

15.0. 

Y 

7 

0.00 

0.05 

0.  10 

0.15 

0.20 

0.25 

0.30 

0.35 

4 

0.000 

0.00 

0.00 

398.18 

585.08 

595.49 

593.58 

613.42 

548.17 

0.050 

0.00 

278.30 

460.38 

575.67 

582.87 

616.81 

604.09 

535.07 

0. 100 

0.00 

376.66 

492. 84 

571.37 

546.13 

573.83 

562.43 

483.78 

0.150 

503.29 

524.27 

532.83 

521.46 

532.11 

517.93 

516.96 

380.49 

0.200 

554.37 

572.62 

564.35 

587.42 

594.77 

536.49 

459.71 

303.94 

0.250 

587.41 

575.34 

591.92 

582.61 

548.57 

397.03 

0.00 

0.00 

0.300 

540.67 

571.90 

570.99 

527.95 

451.21 

311.37 

0.00 

0.00 

0.350 

424.34 

416.06 

416.06 

378.04 

255.07 

0.00 

0.00 

393.54 

w(  I  * J  «  K) 

FOR  SWIRL  VANE 

BLAOE  . 

MGLE  15.0. 

Y 

7 

0.00 

0.05 

0.  10 

0.1  5 

0.20 

0.25 

0.30 

0.35 

L 

0.000 

0.00 

0.00 

255.61 

279.70 

259.38 

3.40 

4.  17 

3.40 

0.050 

0.00 

142.03 

245.  59 

309.46 

268.48 

235.71 

252.01 

87.49 
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0.100 

0.00 

239.  76 

261.95 

302.33 

273.77 

265.67 

262.45 

133.91 

0.150 

251.97 

244.74 

270.28 

289.79 

281.93 

290.19 

262.22 

163.39 

0.200 

270.63 

271.36 

282. 59 

284.51 

269.90 

259.83 

252.30 

220.99 

0.2  50 

286.68 

25  5.36 

285.29 

270.14 

281  .43 

240.79 

0.00 

0.00 

0.300 

268.72 

261.91 

268.  11 

217.91 

191.69 

208.53 

0.00 

0.00 

0.  3  50 

189.84 

185.11 

182.47 

238.24 

181.19 

0.00 

0.00 

0.00 

RY(  It j 

»K )  FOR  SWIRL  VANE  3LA0E 

ANGLE 

15.0. 

Y  0.00  0 

.05  0. 

10  0.15 

i  0.20 

0.25 

0.30  0 

.35 

2 


0.000 

0.000 

0.000 

0.013 

0.023 

0.025 

0.021 

0.027 

0.021 

0.050 

0.000 

0.000 

0.009 

0.027 

0.022 

0.028 

0.027 

0.023 

G.  100 

0.000 

0.00  0 

0.015 

0.025 

0.023 

0.028 

0.027 

0.022 

0.150 

0.018 

0.021 

0.024 

0.023 

0.027 

0.023 

0.025 

0.018 

0.200 

0.020 

0.022 

0.021 

0.025 

0.028 

0.028 

0.024 

0.011 

C.  250 

0.024 

0.021 

0.025 

0.025 

0.024 

0.024 

0.000 

0.000 

0.300 

0.028 

0.02  7 

0.02? 

0.026 

0.024 

0.015 

0.000 

0.000 

0.350 

0.02  5 

0.024 

0.024 

0.021 

0.010 

0.000 

0.000 

0.000 

RYdfJvKI  FOR  SWIRL  VANE  BL4DE  ANGLE  15.0. 

Y  0.00  0.05  0.10  0.15  0.20  0.25  0.30  0.350 

Z 

0.000  0.  0.  5.  5.  4.  0.  0.  0. 


0.050 

0. 

20. 

15. 

8. 

3. 

1. 

2. 

0. 

0.100 

0. 

24. 

14. 

8. 

4. 

2. 

2. 

0. 

0.150 

14. 

12. 

10. 

9. 

5. 

4. 

5. 

0. 

0.200 

13. 

12. 

12. 

9. 

5. 

3. 

2. 

3. 

0.250 

13. 

10. 

10. 

8. 

7. 

0. 

0. 

0. 

0.300 

13. 

11. 

11. 

7. 

4. 

0. 

0. 

0. 

0.350 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

RZU.JtK)  FOR  SWIRLVANE 

BLADE 

ANGLE 

15.0. 

7 

Y  0.00 

0.05 

0.10 

0.  15 

0.  20 

0.25 

0.30 

0.35 

4 

0.000 

0.000 

0.000 

0.012 

0.028 

0.027 

0.030 

0.028 

0.022 

0.050 

0.000 

0.011 

0.023 

0.023 

0.028 

0.027 

0.026 

0.013 

0.100 

0.000 

0.021 

0.021 

0.024 

0.  021 

0.020 

0.019 

0.012 

0.150 

0.020 

0.019 

0.018 

0.018 

0.015 

0.012 

0.014 

0.004 

0.200 

0.025 

0.026 

0.026 

0.026 

0.024 

0.014 

0.00  7 

0.004 

0.250 

0.027 

0.027 

0.027 

0. 02  5 

0.  021 

0.000 

0.000 

0.000 

0.300 

0.016 

0.021 

0.021 

0.014 

0.  005 

0.000 

0.000 

0.000 

0.3  50 

0.000 

0.0  00 

0.000 

0.000 

0.  000 

0.000 

0.000 

0.021 

ALFA!  I » J »  K ) 

FOR  SWIRL  VANE  BLADE  ANGLE  15 

.0. 

Y  0 . 00  3 

.05 

o.io  a 

1.15  0.20  < 

3.25  ( 

3.30  0 

1.35 

I 
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v_r  *.*  >  *.* 


0.000 

0. 

0. 

42. 

27. 

21. 

20. 

20. 

18. 

0.050 

0. 

0. 

34. 

25. 

25. 
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